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WHAT  DO  FLUCTUATION  MEASUREMENTS  SAY  ABOUT  D = Y/k2 
AND  SHEAR  STABILIZATION  OF  DRIFT  WAVES? 

I.  Fluctuation  Measurements  Concerning  D = Y/k2 
A.  Introduction 

The  well-known  estimate  of  the  turbulent  diffusion  coefficient 

\ 

D < “T  , CD 

k2 

x 

(where  Yfc  is  the  growth  rate  and  k is  the  wave  number  of  the  unstable 
waves.)  can  be  derived  in  a variety  of  ways.  One  way  is  to  assume  that 
the  waves  cause  diffusion,  add  a term  - D^n  to  the  ion  continuity 
equation  and  adjust  D so  as  to  produce  marginal  stability.  The  theories 
of  Dupree1,  Weinstock2  and  others3  might  be  regarded  as  more  sophisti- 
cated versions  of  this  procedure. 

However,  the  most  popular  derivation  of  Eq.  (1)  is  that  of 
Kadomtsev4.  This  begins  by  assuming  that  non-linear  effects  saturate 
the  growth  of  the  instability  when 


k2  <n2>  $ 
x 


dn 
o 

dx 


(2) 


(Simple  slab  geometry  is  assumed  with  macroscopic  quantities  varying  in 


the  x direction,  n is  the  mean  and  n the  fluctuating  ion  density.) 
o 

Note:  Manuscript  submitted  November  20,  1978. 
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From  a quasi-linear  treatment  of  the  ion  continuity  equation  one  obtains 
for  the  particle  flux 


<s 


(3) 


If  the  estimate  (2)  for  (rt2)  is  now  used  in  Eq.  (5)  (notwithstanding 

the  inconsistency  of  this  procedure)  Eq.  (1)  follows.  This  estimate  is 

the  basis  of  many  numerical  transport  calculations5. 

Clearly,  the  above  argument  is  not  entirely  convincing  and 

counter-arguments  have  been  given6.  Also  Dupree  and  TetreaultT  have 

recently  reexamined  the  basis  of  the  strong  turbulence  theory.  Their 

conclusion  is  that  D ~ v/k2  can  only  be  valid  if  the  plasma  is  very 

weakly  unstable.  (Their  condition  for  stabilization  by  this  mechanism 

is  v/u)  < (k  v./u))2,  where  k is  the  wave  number  parallel  to  B and  v.  is 
ii  i * ii  ~ i 

the  ion  thermal  velocity.)  Nevertheless  it  is  widely  believed  that 
there  is  a comparatively  simple  relationship  between  diffusion  and 
fluctuation  levels,  that  Eq.  (1)  gives  an  upper  limit  to  the  turbulent 
diffusion,  and  that  Eq.  (2)  gives  an  upper  limit  to  the  level  of  density 
fluctuations.  Evidence  can  be  adduced  in  support  of  these  beliefs8’9. 

In  the  next  section  we  show  that  there  is  also  evidence  against  them. 

B.  Evidence 

It  is  convenient  to  rewrite  Eq.  (2)  in  the  form 

— * (k  L )*1  , (10 

n x n * 

o 

where  is  the  density  scale  length.  We  now  discuss  some  experimental 
evidence. 
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(a)  Okabayashi  and  Arunasalam10  report  studies  of  drift-wave 
turbulence  in  the  FM-1  spherator.  In  their  low  shear  case,  with  the 
plasma  in  a state  of  strong  turbulence,  they  found  the  fluctuation  level 


to  be 


~ 0.15. 


For  (k  L )-1  ( (k.a)-1  in  their  notation)  they  found 

x n v 

(k  L )_1  ~ 0.02  to  0.05,  (6; 

' x n 

so  the  fluctuation  amplitude  is  five  to  seven  times  larger  than  the 
estimate  of  Eq.  (4).  This  translates  (using  Kadomtsev' s prescription) 
into  a diffusion  coefficient  of  between  25  and  50  times  larger  than 
given  in  Eq.  (1). 

(b)  Navratil  and  Post11  observed  large  amplitude  density 
fluctuations  in  a levitated  octupole  which  they  interpret  as  due  to 
drift  waves.  For  the  parameters  of  this  experiment  (Helium  plasma, 

B = 750G,  = 0.2  ev,  k^p^  ~ 0.5)  or  from  the  quoted  wavelength  (5  cm) 


we  infer 


kx  — 2.5  cm"1  . 


The  authors  do  not  state  the  density  scale  length  but  from  their  Fig.  2 
and  previous  descriptions  of  the  device12*13  a reasonable  inference  is 

L ~ 5 cm.  (8] 

n 


(k  l r1  ~ 0.08, 

x n 


whereas  the  measured  fluctuation  level  was 


— ~0.2.  (10) 
n 
o 

The  discrepancy  is  not  as  significant  in  this  case,  but  the  remarkable 
result  from  this  experiment  is  that,  notwithstanding  the  high  fluctua- 
tion  level,  no  anomalous  diffusion  was  observed!  (More  precisely,  the 
transport  estimated  on  the  basis  of  quasi-llnear  theory  and  the  observed 
fluctuation  level  is  — 103  times  greater  than  the  observed  upper  bound 
to  the  anomalous  transport. ) 

(c)  Saison,  Wirnmel  and  Sardei14  solved  numerically  in  two 
dimensions  the  non-linear  trapped- fluid  equations  used  by  Kadomtsev 
and  Pogutse15  as  a description  of  the  trapped-ion  mode.  They  compared 
their  numerically  obtained  values  of  the  diffusion  coefficient  (at  late 
times,  when  the  instability  saturates)  with  Djrp  » the  diffusion  coeffi- 
cient suggested  by  Kadomtsev  and  Pogutse  for  this  instability.  (DKp 
is  in  fact  equal  to  v^/ 2k2.  Linear  theory  does  not  provide  a k^. 
Kadomtsev  and  Pogutse15  assume  approximately  isotropic  turbulence, 
k2  w k2  k2/2,  and  are  followed  in  this  assumption  by  Duchs  et  al.5.) 

Although  Saison  et  al.  found  that  the  density  fluctuations  in 
the  steady  state  were  an  'ordered  pattern'  rather  than  'turbulence',  the 
calculated  diffusion  coefficients  were  up  to  20  times  higher  than  D^p 
and  the  scaling  was  Bohm-like. 

II.  Shear  Stabilization  of  Drift  Waves 
A.  Introduction 

The  basic  theory  of  shear  stabilization  of  drift  waves  is 

le-is.  ^ magnetic  field  is  in  the  z direction 

1+ 


simple  and  physical 


■ 
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and  the  density  gradient  is  in  the  x direction,  the  simplest  theory  of 

drift  waves  gives,  of  course,  cu  = k v_  = k cT  /eBL  (n  L _1  3 dn  /dx), 

yD  ye  non  o ’ 

where  the  notation  is  standard.  If  one  allows  for  ion  inertia  both 
perpendicular  and  parallel  to  the  magnetic  field,  the  dispersion  rela- 
tion becomes 


, a 2 k2^  / k v. 

- k P + z s + [ _j_] 

ur  \ <u 


D - 1 


k2^2  \ _ Q 

y i / 


(11) 


where  p2  = Tg/M*^  s v^/cu2^.  an  Inho^Seneous  system  the  dispersion 
relation,  Eq.  (11),  becomes  a second  order  differential  equation  for  the 


perturbed  potential.  If  x = 0 is  the  point  of  maximum  density  gradient, 
d2  x2 

then  k2  = jj-g-  , k2  = k2  — (where  Lg  is  the  shear  length)  and  vQ  = 
x2  7 s 

vDQ  (1  - ^g-),  say.  (Typically,  ^n~Ln*)  The  equation  for  cp  then  becomes 
n 


K & + 


{A  • a 
■ 1 • » = 


(12) 


assuming  ou  « kyvDo*  Also  we  have  assumed  that  some  effect  (i.e.  inverse 
Landau  damping)  gives  rise  to  a growth  rate  V. 

Equation  (12)  displays  two  fundamentally  different  types  of 

/ v2  \ 

behavior  depending  on  the  sign  of  ( -g-  S - r?  ) . If  it  is  negative, 

\ W's  n / 


localized  modes  can  exist17  and  the  dispersion  relation  is 
U)  = iV  + 


k v 
V DO 


1 + k2P2 


(13) 
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where  n is  an  arbitrary  integer.  Thus  as  long  as 


Finally,  even  if  there  are  no  unstable  eigenfunctions,  it  has 
been  shown  that  drift  waves  can  be  convectively  unstable  and  that  in- 
creasing the  shear  can  stabilize  these  convective  instabilities  by  re- 
ducing the  number  of  growth  lengths  in  the  unstable  region18. 
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Recent  theory  has  diverged  from  this  simple  picture  in  two 
different  ways. 

(a)  In  toroidal  geometry,  variations  in  curvature  drift  terms 
may  lead  to  substantial  reduction  or  even  complete  elimination  of  shear 
damping19*20  (Divergence  A). 


(b)  In  slab  (or  cylindrical)  geometry,  and  in  the  absence  of 


current  along  the  field,21*22  the  electron  response  in  the 

L 

(u/k  > v renders  all  drift  waves  stable  as  long  as  — > 
II  e L„ 


irtw  region 


B.  The  Evidence 

There  are  now  at  least  two  experiments10*23  in  toroidal 

geometry  where  divergence  A might  be  expected  to  apply  but  which  show 

shear  stabilization  of  drift  waves.  In  ref.  23  the  fluctuations,  which 

were  identified  as  drift  waves,  were  observed  to  be  strongly  excited, 

with  amplitude  independent  of  shear  length,  for  Lg  > X^  Vg/VD*  As 

soon  as  L decreased  to  the  point  where  the  above  inequality  reverses , 
s 

there  occurred  a steady  decrease  in  fluctuation  amplitude  with  increas- 
ing shear. 

Secondly,  in  the  experiment  of  Okabayashi  and  Arunasalam10, 

generally  L /L  « p./X  . These  authors  found  that  (i)  increasing 
ns  in 

the  shear  led  to  a decrease  in  the  fluctuations,  (ii)  the  experimentally 
measured  ion -mass  dependence,  and  also  the  absolute  magnitude,  of  the 
shear  marginal  stability  condition  was  in  reasonable  agreement  with  the 
predictions  of  Pearlstein  and  Berk16,  and  Rutherford  and  Frieman18,  and 
(iii)  with  high  and  moderate  shear  the  unstable  waves  could  be  readily 
identified  as  drift  waves. 
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In  slab  geometry  there  is  at  least  one  experiment24  demon- 
strating shear  stabilization.  However  here  L /L  ~ p./X  , so  it  is 

n s I n 

difficult  to  know  whether  stabilization  is  due  to  the  fact  that  elec- 
trons very  near  to  the  point  at  which  k = 0 stabilize  the  modes,  or 
whether  the  more  conventional  theories  are  operative. 

However,  in  the  particle  simulations  of  Lee,  Kuo  and  Okuda25 
the  density  profile  is  initially  exponential  and  therefore  Xr  = °°.  The 
geometry  Is  plane  slab  so  that  divergence  B might  be  expected  to  apply. 
Nevertheless  these  authors  definitely  see  growth  of  unstable  drift 
waves  even  at  fairly  large  values  of  shear.  (Furthermore  they  have 
verified  that  there  is  no  wave  reflection,  either  from  internal 
turning  points  or  from  external  walls.)  Also  they  state  that  the  shear 
marginal  stability  condition  is  reasonably  well  approximated  by  the 
standard  theory  of  Perlstein  and  Berk16. 

III.  Conclusions 

On  the  issue  of  nonlinear  stabilization,  recent  work  shows  that 
fluctuation  levels  can  be  significantly  higher  than  that  suggested  by 
Kadomtsev4  and  diffusion  can  be  significantly  larger  than  y/k2.  Also 
there  does  not  seem  to  be  any  simple  general  relation  between  diffusion 
and  fluctuation  levels. 

On  the  other  hand  recent  results  show  that,  in  spite  of  theoretical 
uncertainties,  the  simplest  theory  of  shear  stabilization  works  fairly 
well  in  both  slab  and  toroidal  geometry. 
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